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1 Introduction 



In the summer of 1995, Dominique Vautherin came to Kyoto and stayed three months 
Q_i! at the Yukawa Institute of Kyoto University where he delivered a series of very infor- 
mal, pedagogical lectures on the application of variational methods to quantum field 
theories, jj]] This initiated our long lasting enjoyable collaboration on the subject 
and Yasuhiko Tsue joined the force later. Throughout this collaboration Vautherin 
was always enthusiastic to explore new physics problems and always came up with 
new innovative ideas, much of which I suspect have origins in his expertise in nuclear 
many-body theory. We are very sorry that this fruitful collaboration with Vautherin 
came to end so soon by his untimely death and that we are no longer able to be inspired 
by his deep physics insights, charmed by his elegance, and most of all cheered up by 
his very warm presence. I like to dedicate these lectures to the memory of Dominique 
Vautherin from whom I learned most of the material presented below. 

First, I like to say a few words about physics motivations of our works. In recent 
years, several authors have constructed non-trivial time-dependent solutions of classical 
field equations of effective meson fields. |2|, [§ Such solutions are relevant in considering 
the fate of defects which might be produced in dynamical order-disorder phase tran- 
sitions in ultrarelativistic nucleus- nucleus collisions [|J or in the evolution of very early 
universe. H We like to study the effect of quantum fluctuations which has been ignored 
in these classical analyses. We may consider the classical solution as a collective mode 
of underlying microscopic degrees of freedom. It is well-known in nuclear many-body 



* Lectures delivered at the NATO Advanced Study Institute on QCD Perspectives on Hot and 
Dense Matter, Cargese, Aug. 6 - 18, 2001. 



1 



theory that the internal microscopic motions of nucleons are influenced by the motion 
of the nuclear mean field such as rotation or vibration. One may expect similar effects 
would arise in the field theory. 

Our method is based on the Schrodinger wave functional representation of the 
quantum field theories in which the time-evolution of the wave function is explicitly 
considered. Another advantage of using the Schrodinger picture is that one can intro- 
duce the variational method to construct approximate but non-perturbative solutions 
of the problem as in similar problems in quantum mechanics. The extension of the vari- 
ational principle for time-dependent wave function has also been developed. || 0] The 
method can also be extended to deal with statistical ensembles described by Gaussian 
form of density matrices. || 



2 A Simple Exercise in Quantum Mechanics 

We start with a very elementary example in quantum mechanics which all students 
learn in an introductory course of quantum mechanics: the harmonic oscillator problem. 
Although our interests lies in the time-dependent solution of quantum field theories, it 
would be instructive to first study the time-dependent solution of this simple exactly 
soluble problem in order to gain some physical insights into our more difficult problems 
since quantum field theories are nothing but an assemble of infinite number of coupled 
harmonic oscillators. 

Let us consider a simple one-dimensional harmonic oscillator whose Hamiltonian is 
given by 



( 2 i 2 2\ 

-{p + u x J 



# = -(y/ + ^./-). <■''.: 

where we set the mass of the particle m = 1 for simplicity. We look for solutions of 
the time-dependent Schrodinger equation: 

z JU(i) = H *(t) 

with the quantization condition: [x,p] = i. 

There are two ways to solve this problem. The easier one is the algebraic method 
in which one introduces "creation" and "annihilation" operators: 

a) = —== (— ip + cox) , a — , — (ip + cox) 



which satisfy a usual commutation relation: a,a^ =1 and diagonalize the Hamilto- 
nian: H = co (a^a + |J . The ground state of the Hamiltonian is given by the condition: 
a\0) = 0. This method is usually transcribed to the quantization of fields and one ob- 
tains particle creation and annihilation operators as basic building blocks in describing 
physical processes. In this approach, one usually does not refer to the wave function 
of the system explicitly but instead concentrates on amplitudes of particular process 
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for given initial and final states specified by the particle number and other quantum 
numbers. 

Alternatively, one can solve the problem by the analytic method in which one 
expresses the Schrodinger equation in coordinate representation with the differential 
operator p = —id/dx and solve the resultant second order differential equation. The 
eigenfunction of this differential equation are Hermite's polynomials with the Gaussian 
ground state wave function: 

/uA 1 / 4 i 2 

tfofc) = (x\0) = (-) e-^ 2 (2) 

It is straightforward to write down the ground state wave functional of the quantized 
free field <p(x). For one component scalar field theory, it becomes just a product of the 
Gaussian wave functions of normal modes each specified by the momentum with the 
oscillator frequency &/. = \/m 2 + k 2 : 



V [tp(x)] = Afexp 



Z k 



with a proper normalization condition. 

Time-dependent variational wave function: 

Now let us consider time-dependent solutions of the harmonic oscillator. We first 
modify the ground state wave function by adding extra complex phase factor ^(po-iuxojx _ 
One then obtains the Gaussian wave function with its the center shifted 



ty(x,t;x ,p ) = exp 



1 2 

-~u(x-xo(t)) +ipo(t)x 



(3) 



Here po(t), Xo(t) are time-dependent parameters which are to be determined by impos- 
ing that the above function is a solution of the time-dependent Schrodinger equation. 
Here we use the variational method to derive the equations of motion of po(t), xo(t). 

The Schrodinger equation can be obtained by imposing a stationary condition 5S = 
for the action: 

sm = fdt(*(t)\H-i-\*{t)) 

with respect to variation of the wave function (ty(t)\. For the variational wave function 
in the form of |^, the integrand can be computed easily with (\I/|p 2 |\l/) = pjj, (\E f |a; 2 |\I / ) = 
xl |*) = -p x : 

d 

(^\H-i— 1$) =pI + uj 2 xI+pqx . (4) 

Taking stationary conditions of the action with respect to two time-dependent param- 
eters, 5S/5po = 5S/5xo = 0, one finds 

po = -uj 2 x , x = po (5) 
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which are just the classical equation of motion of the harmonic oscillator. The center- 
shifted Gaussian wave function may be obtained, beside the phase factor e~ ip ° x , by 
operating a translation operator on the ground state: 



\0;x ,p ) ~ exp[-ix p]\0) = exp 



|0>. 



//A 1 
* sq .0M) = ( - ) exp 



It thus describes a phase-coherent mixture of infinite number of excited states \n). The 
corresponding states in quantum field theory are called coherent states. One thus sees 
that classical field equations arises if one describes the wave functional of the quantized 
fields in terms of coherent states. 

Keeping Gaussian form of variational wave function one can go one step further by 
introducing additional phase factor quadratic in x. The wave function now takes the 
form of 

1/4 r I 

■-(// + ia)(x - x ) 2 - ip x (6) 

where we have introduced two new real parameters, /i, a, characterizing the width of the 
wave function, which we shall consider as time-dependent variational parameters. The 
time-dependent parameter /z(t) is related to the quantum fluctuation of the position of 
the particle around its mean x = (x) by ((x — x ) 2 ) = l//i(t). The another parameter 
a is related to the rate of change of /i as we shall see below. This modification of 
the wave function thus describes the breathing motion (squeezing and stretching) of 
the wave function centered at x (t). In the quantum field theory, these generalized 
coherent states are called squeezed states. 

With the new variational wave function the action integrand becomes 

WW - ^ t Mt)U. = \{fo + " 2 4) + \{^ + 

Taking the stationary conditions with respect to the variations of /i(t), cr(t), in addition 
to classical parameters, xo(t),po(t), one finds 

jx = 2a fi, a = a 2 + uj 2 — fi 2 , (7) 

while the equations of motion of xo(t),po(t) are unchanged. We see that the imaginary 
part a of the Gaussian width parameter plays a role similar to the velocity of the 
motion of the center of the Gaussian x§. For the harmonic oscillator Hamiltonian, 
coherent states and squeezed states are exact solutions of the Schrodinger equations. 
The classical motion and the quantum fluctuation decouple in this exactly soluble 
problem. This is not the case, however, when the potential is not quadratic in x. 

Anharmonicity: 

To illustrate the effect of non-harmonic part of the potential, we add a term quartic 
in x in our Hamiltonian: 
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In the scalar field theory, this corresponds to adding a (p 4 self-interaction term. This 
anharmonic term generates new terms in the integrand of the action: 



(V(t)\Hr\*(t)U. = - 




which cause coupling between the classical motion of the mean of particle position and 
the quantum fluctuation around it: 

xo = Po, Po = -u> 2 x - ^xl - ^-x , 

6 4/t 

2 9 2 ^ ^ 2 

H = 2cr//, a = a + uj — n + h -x 

4/i 2 

Let us see how the time independent solution (the ground state) is modified by the 
anharmonic term. The conditions Xq = /io = demand that po = a = 0, while the 
remaining two conditions po = cr = determine the values of xo and fi. There are 
two types of solutions depending on the sign of uj 2 . When uj 2 > we have a "normal" 
solution centered at the origin xo = but the width is modified slightly as determined 
by: 

/i 2 = uj 2 + — 
4/i 

This equation may be called the gap equation for the reason to be discussed in the next 
section. If uj 2 < 0, the potential has double minima at rc min . = ±^6/~\uj and the x = 
point becomes local maximum of the potential. In this case, there are two "symmetry 
breaking" solutions centered at the two solutions of 

where the value of the width parameter /t is determined self-consistently with the 
modified gap equation: 

We note that the position of the center of our Gaussian variational wave function shift 
slightly toward the origin from the position of the minimum of the potential. 



3 Scalar Field Theory: tp A model 

Having been warmed up by a much simpler problem in quantum mechanics, it is now 
our task to transcribe the result to a problem in quantum field theories. We consider 
first a prototype scalar field theory with <p 4 self-interaction. 
The Hamiltonian density of the theory is given by 

W(x) = ^ 2 (x) + \ (Vv?(x)) 2 + ^ 2 (x) + ^ 4 (x) , (8) 
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Table 1: Correspondence between quantum mechanics and scalar field theory 



quantum mechanics 


scalar field theory 


X 


p(x) 


cu 2 


ml — A ( or ml + k 2 ) 


4 X 

x (t),p (t) 




V?(x,t),7r(x,t) 




|G- 1 (x,y,t),2S(x,y,t) 




G~ 1/2 £G~ 1/2 



where 7r(x) is an operator conjugate to the field <^(x) and is expressed by a func- 
tional derivative —i8/5(p(x). [] We write a Gaussian time- dependent variational wave 
functional for this Hamiltonian formally as 

* b(x)] = Af exp fi<7f|p - £>} - - <p\jg + iE\<P - , (9) 

where G, E, 7f define respectively the real and imaginary part of the kernel of the 
Gaussian width and its average position and momentum. We have used the short 
hand notation (tv\(p) = J 7r(x, t)(p(x)dx. Although it looks a little horrible, it is just 
a straightforward generalization of The correspondences between the previous 
quantum mechanical example and the present case are summarized in the Table 1. 
The equations of motion are found to be 

if = -7T, 

7T = (-A + m 2 + ^ 2 (x) + fG(x,x))^, 

(10) 

G = 2(GE + EG), 

t = lG- 2 -2S 2 -i(-A + m 2 + |^ 2 + |G(x,x)). 

where it is understood that (p and 7f denote vectors with components <^(x) and vf(x), 
G and E denote matrices with matrix elements G(x, y) and E(x, y), and the matrix 
product GE is given by / dzG(x, z)E(z, y). 

For the vacuum, we should have time-independent solution so that tt — E = 0, 
y?(x) = ipo, and 

g(x,y. t )=G (x-y) = / rf k i g^. (11) 

1 A care must be taken to give precise mathematical meanings for these expressions ||, but we do 
not go into such problems here. 
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where the effective mass // is determined self-consistently by the non-linear integral 
equation 

/i 2 = ml + ^Go(O) + ^ . (12) 

which is usually called the gap equation because it determines the mass gap self- 
consistently. In the symmetric phase the expectation value of the field (fo vanishes 
while in the symmetry broken phase it must be such that 

™S + ^G o (0) + ^g = 0. (13) 

This last equation implies that /i 2 =A<£>q/3. The equations of motion (|II]) determine 
the time-evolution of the variational wave functional (|9|). 

We comment on two well-known problems which are absent in quantum mechanics 
but are characteristic in quantum field theory: divergences and covariance. In the above 
expression the momentum integral in Gq(0) is quadratically divergent. As well-known, 
this divergence originates from couplings of infinite degrees of freedom (all momentum 
modes) in the interaction of local fields. In the renormalizable field theories, these 
divergences (or cut-off dependence) may be absorbed into an appropriate redefinition of 
finite number of physical parameters, such as the mass and the coupling constant. This 
renormalization procedure works at least order by order in the power series expansion 
in terms of the coupling constant. In our non-perturbative calculation scheme with 
the variational method, we may absorb the cut-off dependence into the mass and 



the coupling constant. [[10], |TT[ We encounter, however, well-known triviality problem: 
renormalized coupling constant A# becomes zero as one sends the cut-off A to infinity 
keeping the original coupling A positive for the stability of the ground state. Despite 
this well-known pathology, the model still makes sense as an effective theory with the 
finite cut-off. 

The another problem is the lack of manifest covariance in our formulation: the 
time coordinate t have been treated differently from spatial coordinates throughout 
calculations. This is an old problem which was originally solved in QED by Tomonaga, 
Schwinger, Feynman and Dyson who developed the covariant perturbation theory based 
on the interaction representation. Vautherin invented a new ingenious trick to derive 



manifest covariant form of the equations of motion which I shall now describe. |12, 14}] 

Mean field equations in the Hartree-Bogoliubov form: 

A key ingredient of his method is the reduced density matrix defined by 

V ; V M x My)> -«<7r(x)p(y)) - 1/2 J v ; 

where = (p — ip , n = n — 7f, and expectation values are calculated with the Gaussian 
functional ^(t) and is given by 

KA={ - 2iG ^ G \ n ^ 
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Using (|TOD, one can show that the equation of motion of the reduced density matrix 
can be cast into the Liouville-von Neumann form 

iM = [H,M], (16) 
where the generalized Hamiltonian 7i is given by 

K=(?J). (IT) 

with 

r = -A + m^ + ^ 2 + ^G(x,x). (18) 

This form of equations is known in many-body theory as the time-dependent Hartree- 
Bogoliubov equations. 

The reduced density matrix satisfies Ai 2 = |l so that it has two eigenvalues of 
±1/2. We write the n-th eigenvector of M. with eigenvalue 1/2 as (w n (x), t> n (x)): 

m {::)-\(::\ ^ 

then one can show that vectors (u*(x), — t>*(x)) give eigenvectors for eigenvalue —1/2. 
The u and v components of eigenvectors are called mode functions. With these eigen- 
vectors the reduced density matrix has a spectral decomposition as 



Z n>0 



v n , u n ) +| "I (-«„, u r , 



(20) 



The Liouville-von Neumann equation (p~6|) can be rewritten in terms of the mode 
functions as 

*(;)-(? :)(:;)■ <*> 

Eliminating v n and inserting flTJD we obtain a modified Klein-Gordon type equation 
for the mode functions 

(n + m 2 + ^ 2 + ^G(x, x)^ u n = 0, (22) 
where the spectral representation ( PUD implies 

G(x,x) = (x|G(t)|x) = i5:K,(x,t)| 2 . 

^ n 

To write the equation of motion fully covariant way, we introduce a Feynman propa- 
gator in terms of the mode functions 

(x\S\y) = 0(x -y )J2n>oK(^ x o) u n(y,yo) 

(23) 

+ d(yo - xo)En<oK(^ x o)un(y,yo)■ 
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so that 



{x\G(x )\'x) = (x\S\x) 



(24) 



Then we finally arrive at a fully covariant, self-consistent equations of motion for the 
mode functions |]Ti2"|, II | which reproduces the same equations obtained earlier by the 
functional integral method. [|Tj| 

This way of writing the equations of motion also paves a way to generalize the 
calculation at finite temperatures. 



4 Statistical Ensembles 

Foregoing discussions are limited to evolution of a single coherent (Gaussian) state. In 
realistic physical situations, we are more interested in the evolution of the statistical 
ensemble which is described by the density matrix: 

p(t) = J2\Mt))pn(t)(M-t)\ (25) 

n 

where p n (t) is a probability distribution specifying the ensemble hence It should satisfy 
J2nPn(t) = 1. The expectation value of an observable O is given by (O) = TrpO and 
the statistical entropy S is given by 

S = -Trplnp = -J^Pn^Pn (26) 

n 

Hence if p n is time-independent then the entropy is conserved. In equilibrium at 
temperature T = 1/(5, the density matrix is given by the canonical ensemble: p n q '(t) = 
e~ f3En /Z where Z = Y^n e ~ /3En which maximizes the entropy S(t) under the condition 
of fixed expectation value of the energy E = (H). The density matrix with time- 
independent p n obeys the Liouville-von Neumann equation: id t p = [H, p]. 

The variational method has been extended for the time- dependent density matrix 
by Eboli, Jackiw and Pi. || Without going into detail, we illustrate the essence of 
their method in terms of simple quantum mechanical example with harmonic oscillator 
Hamiltonian ([!]). We first introduce the coordinate representation of the density matrix 
by 

p(x,y;t) = (x\p(t)\y) = Y,*n(x,t)Pn(t)K(y,-t) (27) 

n 

and observe that with single Gaussian variational wave function the density matrix is 
just a product of two Gaussian: p(x,y;t) ~ e ~^ x +v >' 2 . One can show that in the 
other extreme limit of thermal equilibrium, the density matrix p(x,y;t) can also be 
expressed by a mixed Gaussian form. This is so because the equilibrium density matrix 
Peq.(P) = e~^ H obeys the imaginary time Schrodinger equation (the Bloch equation): 
—d T p eq .{ T ) = Hp eq . Using p eq (x,y) = (x\U(—ir)\y) and well-known path-integral 
expression of the matrix elements of the unitary evolution operator U(t) = e~ ttH , one 



9 



finds 



Pe q .(x,y) 



ui 



1/2 



exp 



UJ 



2 sinhuft 



|(x 2 + y 2 ) coshcj/5 — 2xyj 



(28) 



which is again Gaussian with an extra term containing a product of two coordinates 
xy. 

For more general mixed states we introduce a generalized Gaussian density matrix || 



p(x,y;t) = A/" exp 



ipo(x — y) — — [x 2 + y 2 — 2£xy) + i—(x 



.a 



y 



(29) 



where we have used the time-dependent shifted coordinates: x = x — xq and y = y — Xq. 
Equations of motion of the time-dependent parameters of the generalized Gaussian 
density matrix (|29"D can be derived from the Liouville equation of the density matrix 
and we obtain a set of equation similar to ( |5UD with a small modification: 

fi = 2afi, a = a 2 + u 2 - (1 - £ 2 )/i 2 . (30) 

The parameter £ is called the mixing parameter which measures the degree of mixture 
of different pure states in the ensemble; it remains constant for an adiabatic evolution 
of the system. In equilibrium, po = xq = a = and other two parameters are given by 
the specific functions of temperature: 



/i e(? . = ojcoth.ujf3, £ e(? . = cosh 1 u)[3 



(31) 



as indicated by the formula (p8|). Extension of the Gaussian density matrix in quantum 
mechanics to that in quantum field theories is straightforward as indicated in the last 
row of the Table 1. 

The reduced density matrices we have introduced in the previous section for a pure 
Gaussian state can be extended for a mixed state immediately by the replacements: 

(£(x)7r(y)) = Tr(p(t)^(x)7r(y)), (y3(x)y3(y)> = Tr(p(^(x)£(y)), etc. (32) 

and one can derive equations of motion of the reduced density matrix similar to the 
pure state case. In the case of equilibrium distribution, this amounts to introduce a 
factor containing the occupation number in the sum over the mode functions. 



5 Rotating chiral condensate in the 0{N) sigma model 

We briefly mention about an application of the above method to the O(N) sigma model 
which is composed of iV-components coupled scalar fields ip n with continuous O(N) 
symmetry. We expect that this global symmetry of the model is broken spontaneously 
at low temperatures, characterized by non vanishing expectation value of one compo- 
nent of fields, say v? = Trptp ^ 0, and the system exhibits an order- disorder phase 
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transition to a state with (pQ = 0. This is what is expected in QCD where the chiral 
symmetry is broken in vacuum and is expected to be restored at finite temperature. 
This chiral phase transition has been studied by an effective theory with pion and 
sigma fields with 0(4) global symmetry. 

We have applied our method to describe a special kind of time-dependent conden- 
sate which rotates in a subspace of the internal symmetry space: 

where r 2 = ( ^ ^ ^ and = (u>, q) is a four vector specifying the direction of 

rotation: for a spatially uniform rotation in time it is a pure time-like vector; while a 
static condensate with oscillation in space, it is a pure space-like vector. The matrix 
U(x) = exp[i(q ■ a;)r 2 ] can be considered as a "gauge transformation" to the local rest 
frame of the rotating condensate. In the rotating frame, the four derivative which 
appear in the equation of motion of the mode functions is "gauge transformed" to 
Ud^JJ^ = d^ — iq^T^ and the effect of the rotation is seen in this frame as an appearance 
of the apparent "centrifugal force". Indeed, our result of phase diagram for such 
dynamical condensates shows that the amplitude of the chiral condensate with uniform 
time-like rotation increases due to the centrifugal force. This effect has been observed 
in the classical solutions of Anselm and Ryskin ||; our quantum generalization of their 
solution shows that this effect is amplified by the coupling of quantum fluctuations to 
rotations, which the static condensate with spatial oscillations are more suppressed by 



the quantum fluctuations. [H| A phase diagram of rotating condensate was obtained 
in [14| and the damping of the rotating condensate due to the symmetry breaking 
perturbation was computed by the method of the response function in 



6 Outlook 

Vautherin started to work on the variational approach to quantum field theories many 
years ago with Arthur Kerman. They developed many important ideas in their un- 
published works and tried to solve QCD non-perturbatively with their method with a 
hope to gain new insights in the quark confinement problem. [|H|] The Gaussian Ansatz 
for the variational wave functional however has difficulty of breaking the local gauge 
invariance and the projection to color singlet state destroys a nice feature of the Gaus- 



sian wave functional. [17|] Vautherin continued to work on the problem with his students 
and brought a new insight into the problem again introducing a technique developed 
in nuclear many-body theory in his last paper. |Tj| His efforts in this direction may 



be carried over to study the dynamical evolution of the quark-gluon plasma in ul- 
trarelativistic nucleus-nucleus collisions. Vautherin was also interested in the recent 
experimental breakthrough of creating weakly interacting Bose-Einstein condensates 
in well-controlled laboratory environments. Our method can be also applied to such 



11 



problem to investigate the effect of quantum fluctuations which are usually ignored in 
theoretical descriptions. 

I am much indebted to Yasuhiko Tsue as well as to Dominique Vautherin for our 
works quoted above. I thank them for sharing the joy of physics in our collaboration. 
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